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ABSTRACT 

Asymptotic  equations  describing  slowly  varying 
planetary-scale  motions  of  the  atmosphere  over  the 
rotating  earth  are  derived  and  compared  to  the 
geostrophlc  equations  of  motion  In  the  tangent  plane 
approximation.   Since  the  same  approximation  procedure 
Is  used  In  both  the  spherical  and  plane  geometry,  a 
direct  term  by  term  comparison  Is  possible.   The  plane 
approximation  yields  familiar  prognostic  equations 
(initial  value  problem);  but  In  the  more  constrained 
spherical  geometry,  a  more  restricted  forcing-type 
problem  results  If  the  zonal  motion  Is  t Ime - depend en t . 
That  Is,  with  a  given  slowly  varying  or  steady  zonal 
motion.  Is  associated  a  class  of  long-wave  tesseral 
motions.  I.e.  coupled  low-wave-niimber  and  meridional 
motions . 
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1.   Introduction 

In  recent  years  dynamic  meteorologists  have  been 
attempting  to  obtain  approximate  equations  of  motion 
describing  the  large-scale,  nearly-horizontal  flow  of 
the  atmosphere  valid  over  the  entire  rotating  earth. 
These  studies  are  a  continuation  of  the  pioneering  work 
of  C.-G.  Rossby  and  have  been  greatly  stimulated  by  the 
rapid  development  of  large  electronic  computers.   A 
review  article  by  Phillips  [1]  indicates  the  present 
understanding  about  this  problem. 

In  this  paper  we  derive  approximate  equations 
describing  adlabatic  planetary-scale  motions  of  the 
atmosphere  valid  for  asymptotically  large  time  from 
the  viewpoint  of  the  original  fluid  equations .   Although 
this  approximation  procedure  yields  an  ordered  system  of 
linear  partial  differential  equations,  the  coefficients 
are  variable;  and  numerical  calculations  must  be 
carried  out  in  the  future  in  order  to  verify  if  the 
actual  physical  atmosphere  allows  such  motions  to  exist. 

One  of  the  principal  results  which  we  can  expect  to 
obtain  from  numerical  studies  of  these  asymptotic  equations 


What  is  implied  by  these  asymptotic  equations  is  that 
heating  and  dissipation  (not  included  in  the  analysis) 
act  mainly  on  smaller  scale  elements  in  the  atmosphere 
and  do  not  directly  modify  the  large  scale  motions 
except  the  purely  zonal  flow  itself  which  must  have  an 
external  driving  mechanism.   This  is  only  approximately 
true  at  best. 


is  the  determination  of  the  secondary  motions  associated 
with  a  slowly  varying  zonal  motion.   These  secondary 
motions  consist  of  a  low-wave-number  eddy  motion 
closely  coupled  with  a  non-geostrophic  meridional  motion. 
The  novelty  of  these  equations  is  that  they  are  apparently 
valid  over  the  entire  earth,  including  the  poles  and  the 
equator  where  the  usual  geostrophic  approximation  formally 
breaks  down.   However  the  highly  constrained  nature  of 
this  rather  degenerate  system  of  equations  seems  to 
restrict  the  asymptotic  motion  to  slowly-time-varying 
nearly-zonal  flows  combined  with  an  associated  meridional 
motion.   In  §3-1  we  find  that  long-wave  tesseral  motions 
are  eliminated  by  the  approximation  procedure  when  we 
perturb  on  the  atmosphere  at  rest  with  respect  to  the 
rotating  earth.   But  in  §3-2  we  show  how  the  more  complex 
tesseral  motions,  which  are  associated  with  a  prescribed 
slowly  varying  zonal  motion,  can  be  studied. 

The  mathematical  technique  used  here  is  the  same  as 
that  which  yields  the  approximate  first-order  geostrophic 
motion  in  the  tangent-plane  approximation  (Morikawa  [2]). 
The  formal  procedure  is  the  following: 

1)   Scale  the  time  t  with  respect  to  a  small  parameter 
e,  where  e  <<  1;  the  new  time  variable  is  t  =  et. 
The  small  parameter  e  may  be  related  to  the 
strength  or  amplitude  of  disturbances  in  the 


initial  or  boundary  data  or  to  the  frequency  of 
the  slowly  varying  motions. 
2)   Express  all  the  dependent  physical  variables 
as  power  series  expansions  in  e,  e.g.  in  §3.1, 
perturbed  on  the  atmosphere  at  rest  with  respect 
to  the  rotating  earth;   the  zeroth  order  terms 
describe  the  solid  body  rotation  of  a  stable 
gravitating,  compressible  atmosphere, 
stationary  with  respect  to  the  earth. 

The  time  scaling  (or  "time  stretching")  performs  the 
task  of  filtering  out  all  but  the  slowly  varying,  large-scale 
motions  which  are  directly  related  to  the  asymptotically 
large-time  behavior  of  solutions  of  the  original  system 
of  equations . 

Now  we  make  some  conjectures:   A  possible  implication 
of  these  asymptotic  equations  is  that  for  some  large  time 
scale,  the  planetary  meridional  circulation  is  driven  mainly 
by  the  zonal  flow,  but  for  smaller  time  (or  length)  scales 
the  zonal  flow  is  forced  by  the  smaller  scale  motions,  i.e. 
cyclonic  and  anti-cyclonic  motions  which  have  more 
detailed  secondary  flow  structure,  directly  associated 
with  them,  than  the  larger,  cellular  meridional  motions. 
Our  results  here  might  be  regarded  as  being  complementary 
to  both  the  tangent-plane  and  Rossby's  p -plane  approximations. 
That  is,  the  present  asymptotic  analysis  describes 


primarily  some  planetary  motions  in  the  long  wave  limit 
of  zero  wave  number  in  §5-1  and  low-vjave  number  in  §3-2; 
the  plane  approximations  appear  to  have  validity  in  some 
range  of  intermediate  wave  niirribers  between  zero  and 
infinity  (say,  between  five  and  fifty) .   There  is  also 
the  possibility  that  the  actual  atmosphere  allows  a 
related  spaclally  complementary  description  within  the 
framework  of  our  asymptotic  equations.   That  is,  the 
planetary  analysis  may  be  approximately  valid  only  near 
the  equator  and  possibly  near  the  poles;  and  the  tangent 
plane  description  is  approximately  valid  only  in  middle 
latitudes.   A  more  comprehensive  analysis,  which  includes 
a  greater  range  of  atmospheric  phenomena  and  which 

includes  both  approximations  as  special  cases,  is 

* 
needed  to  explain  the  precise  relationship  between  them. 

This  may  require  a  more  general  approximation  procedure 

which  allows  the  introduction  of  at  least  one  additional 

time  (or  length)  scale.   The  direction  in  which  to  search 

for  such  a  formulation  may  lie  in  the  study  of  large-scale 

atmosphic  turbulence;  then  heating  and  dissipation  must 

be  at  least  implicit  in  the  future  analysis,  which  would 

be  a  difficult  task  indeed. 


■X- 

A  thorough  numerical  study  of  the  asymptotic  equations 
described  in  §5.2  may  furnish  some  relevant  information 
about  this  relationship.   However,  the  required  global 
meteorological  data,  preferably  daily  over  a  period  of 
at  least  one  year,  is  apparently  not  yet  available. 


2.   Conservation  of  Potential  Vortlclty  and  the 
Equations  of  Motion. 

In  the  study  of  large-scale  motions  of  the  atmosphere 
In  which  the  rotation  of  the  earth  plays  an  Important  role, 
the  concept  of  potential  vortlclty  advectlon  (or  convection) 
Is  vital.   This  Is  especially  accentuated  In  adlabatlc 
flows  since,  there,  the  potential  vortlclty  Is  conserved 
on  particle  paths.   In  a  frame  of  reference  rotating 
with  constant  angular  velocity  O^  the  potential 
vortlclty  conservation  theorem,  which  Is  attributed   to 
Euler  (1757),  Ertel  (19^2)  and  Rosshy  (19^9),  states 

(1)         1^  {(20  +  I)      (p'Vad  0)}  =  0 

where  d(  )/dt  Is  the  three-dimensional  particle  (or 

->       -> 
material)  derivative  operator,  ^  ^  curl  v  Is  the  relative 

vortlclty,  p  Is  the  density,  and  0  Is  the  potential 

temperature  (directly  related  to  the  entropy  S  =  C  In  0  , 

where  C   Is  the  specific  heat  at  constant  pressure) . 

The  derivation  of  Eq.  (1)  requires  the  explicit 

use  of  the  following  equations  valid  In  the  rotating 

reference  frame: 


In  the  more  recent  literature  an  exposition  of  these 
and  related  matters  Is  given  by  Truesdell  [5]  In  his 
treatise  on  vortlclty;  more  directly  on  meteorological 
subjects,,  we  refer   to  Ellassen  and  Klelnschmldt  [5]. 


1)  Conservation  of  mass  (or  continuity  equation) 

(2)  ll  +  p  div  V  =:  0 

2)  Vorticity  advection  equation  obtained  by  operating 
with  the  curl  (  )  on  the  inoment;Am  equation, 

■,— >         _>      _>  -. 

(3)  -^  +  20xv+p~      grad  p   +  grad  ^  =  0 

where 

dv^    ->  ->   ->    ->        "^    1       i->i2 

(4)  -tt:  =  V,  +  (v  •  grad)v  =  v,  -  v  x  curl  v  +  -p  grad  |v  | 

and  p  is  the  fluid  pressure,  ^  is  the  gravitational 
potential,  and  subscripts  Indicate  partial  differentiation; 
the  resulting  vorticity  advection  equation  is 

(5)  -if  +  (20+  O    dlv  V  -  [  (20+C) -gradlv  +(grad  p  ^)    (grad  p)-0 

3)  Conservation  of  potential  temperature  for  adiabatic 
flow, 

f^^  IT  =  ^ 

where  the  equation  of  state  relates  ©  ,  p  and  p, 

7-1     1 
(7)  ©  =  (7Po^  -P^)  /  P 

where  Pp,  is  the  pressure  at  the  ground,  r  =  a,  and  the  gas 


constant  7  =  C  /C  ,  the  ratio  of  specific  heats. 

For  our  basic  complete  system  of  equations  we  choose 
Eqs .  (1),  (3),  (6)  and  (?)  and  the  appropriate  Initial 
and  boundary  conditions  on  the  surface  of  the  earth  and 
at  Infinity.   A  completely  equivalent  system  of  equations 
Is  given  by  replacing  the  potential  vorticlty  equation  (l) 
by  the  continuity  equation  (2);  however,  the  asymptotic 
behavior  of  the  motion  is  obtained  more  succinctly  by 
suppressing  the  explicit  use  of  Eq .  (2)  (cf.  [2]). 
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3-   Asymptotic  Equations  of  Large-Scale  Planetary  Motions 

The  formal  approximation  procedure  which  we  use  to 
obtain  the  asymptotic  equations  Is  quite  simple.   That  Is, 
the  rules  are  elementary  although  the  complete  physical 
and  mathematical  justification  and  precise  estimates  of 
error  Involved  In  the  approximation  are  not  totally  apparent 
at  the  present  time.   Previously  we  have  used  this  same 
procedure  to  obtain  approximate  equations  describing 
two-dimensional  geostrophlc  vortex  motion  on  a  plane 
tangent  to  the  rotating  earth  (Ref.  [2]).   In  Appendix  I, 
we  derive  the  equivalent  geostrophlc  equations  of  motion 
for  an  adlabatlc,  compressible  atmosphere  over  a  tangent 
plane,  i.e.  In  which  the  horizontal  and  vertical  components 
of  the  earth's  rotation  rate  vector  O  is  assumed  to  be 
constant  and  fixed  to  the  latitude  at  the  point  of 
tangency.   In  Appendix  II  we  show  the  relationship 
between  the  vertical  structure  of  the  atmosphere  (at  rest 
with  respect  to  the  tangent  plane)  and  the  vertical 
variation  of  the  horizontal  length  scale  of  a  singular 
geostrophlc  vortex  line  element.    At  the  present  time 


Apparently  these  equations  are  similar  to  those  first 
derived  by  Obukhov  [5]  and  Monln  [6]  by  a  related  method. 

In  addition,  the  strength  of  the  vortex  is  functionally 
related  to  the  vertical  structure  of  the  static  atmosphere 
and  hence  to  the  horizontal  length  scale. 
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the  motivation  for,  and  validity  of,  the  asymptotic 
equations  must  rest  partly  on  the  successful  derivation 
of  approximate  equations  for  geostrophic-motion  over 
the  tangent  plane,  as  derived  in  these  Appendices.   In 
addition,  further  substantiation  must  await  the  testing 
of  these  asymptotic  equations  by  comparison  with  actual 
atmospheric  motions,  using  numerical  methods  aided  by 
electronic  computers . 

5.1.   Asymptotic  Equations  Obtained  by  Perturbing  on  the 
Atmosphere  at  Rest. 

In  order  to  proceed  with  the  derivation  we  first 
transform  the  basic  equations  (1),  iJ>) ,    and  (6),  which  are 
in  vector  form,  into  equations  in  spherical  polar 
coordinates  and  introduce  the  scaled  time  t  =  et.   Then 
the  independent  variables  are  (0,  (j),  r,  x)  where  0  is 
longitude  increasing  eastward,  ^   is  latitude  increasing 
northward  from  the  equator,  and  r  is  radius  measured 
from  the  earth's  center.   The  velocity  components  corres- 
ponding to  the  coordinates  (0,  ^,    r)  are  (u,  v,  w)  =  v. 
Thus  the  dependent  variables  expressed  as  a  state  vector 
(u,  V,  w,  p,  p,@)  for  the  atmosphere  at  rest  (e  =  0)  is 
(0,  0,  0,  p^   ,  p^   ,  0    )  representing  the  lowest  order 
terms  in  the  power  series  expansion  with  respect  to  e  and 
is  independent  of  time.   The  potential  vorticity  equation 


(1)  is  transformed  Into  spherical  polar  coordinates  using 
the  following  identities:   The  earth's  rotation  vector 
becomes 

(8)  20  =  (0,  2  0cos  (j),  20  sin  (j)) 

and  we  abbreviate  by  calling  the  (f)- component,  20  cos  (j)  =  fp, 
and  the  r-component,  20  sin  (\>  =    f',    the  relative  vorticity 
is 

(9)  C  =  curl  V  =  (e,  T],  0 

=  [^   [w^-(rv)^],  ^  [(ru)^  -  ^3!^]  ,  ^  g^3  ^  [vq-(u  cos  <i))^]}; 
the  gradient  of  the  potential  temperature  is 

(10)  grad  0  =  ( ^  ,  — ^  ,  0  )   ; 


and  the  particle  derivative  operator  on  a  scalar  quantity 
is 

(11)  ^y-  =  £(  )   +  ^^^— r  (  )q  +-  (  )l  +  W(  ) 

^     dt     ^   T   r  cos  9^0   r  ^   (p    ^    'x> 

where  lettered  subscripts  denote  partial  differentiation. 
Also  the  potential  temperature  equation  (6)  is  immediately 
transformed  into  spherical  polar  coordinates  by  Eq.  (11) . 
The  (0,  (j),  r)  components  of  the  momentum  equation  (3)  are 
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^^2)     It   +  ~ r +  pr   cos   ^  +  ^O  (w   cos   ^-y  sm   ^)=0 

(13)     ^  +  ^  + ^ — ^+  ^  +  20  u  sm   (j)  =  0 

p  p  Y>  p 

Mjn      dw        (u    +  V   )     ,   ^r    .    a  g        ^^  i         „ 

^^    '     "dt   ~  ^ r ^       "o~       — 2   "      ^^   '^ ^   4)  =  0 

r  r> 


where  a   Is   the   earth's   radius  and  g   Is   the  gravitational 


-> 


acceleration  at  r  =  a.   For  completeness,  div  v  becomes 

->                    Uq      (v  cos  (|))  I    (r  w)^ 
(15)   dlv  V  =  dlv  (u,  V,  w)  =  r  H T— ^  -f T^ . 

The  next  step  in  the  derivation  Is  to  express  the 
state  vector  (u,v,w,p,p,0  )  as  a  power  series  expansion 
with  respect  to  the  small  parameter,  e, 

(16)  (u,v,w,p,p,0)  =  (o,o,o,p^°^pf°^0f°h  + 

,  (1)  (1)  (1)  (1)  (1)  ^(1). 

+  e^u    ,  V    ,w    ,p    ,p    ,  y    )  + 


^   2,  (2)   (2)   (2)   (2)   (2)  ^(2),  ,^, 
+  e  (u^  '  ,v^    Sw^  Sp^  Sp^   ,0    )+0(i 


The  final  main  step  Is  to  Insert  the  expansion  (I6)  Into 
the  basic  vector  equations  (1),  (3)^  (6)  which  have  been 
transformed  using  Eqs .  (8)  to  (l4)  and  equating  terms 
having  equal  powers  of  e  in  the  resulting  equations  in 
the  usual  way.   Also  the  equation  of  state  (7)  must  be 
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expanded  and  becomes 


(17)  0  -    0^°^  +  £0^^^  +  0(e2) 

7P      P 

^(0)    ,    7-In/-  (0)^l/7  /   (0) 
where   ©^    =  (7  Pq   )  (p    )     /  P    • 

The  first-order  equation  of  state  becomes 

^(1)    ^(1)      fl) 

To  follow  the  prescribed  calculating  routine  outlined 
above  is  quite  elementary  so  we  summarize  the  results  of 
the  approximation  procedure  to  0{e    ): 

1)   Zeroth  order  terms  --  From  the  {Q,^)    components 
of  the  momentum  equation  (12)  and  (13)  we  get  to  lowest 
order  in  e  that  pA   =.  0  and  pi   =  0  respectively;  hence 
p^   (r)  is  a  function  of  r  only,  having  assumed  at  the 
outset  that  p^  '    =  0.      Then  from  the  r-component  of  the 
momentum  equation  (l4)  to  lowest  order  and  the  equation 
of  state  (17),  we  get  p^   (r)  and  consequently  0^"^  (r); 
the  atmosphere  at  rest  with  respect  to  the  rotating 
earth  is  therefore  spherically  symmetric  and  hydrostatic, 

pt°'     2 

(19)  -TOT  +  -H  =  0  • 

p^      r 
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2)  First-order  asymptotic  equations  --  Prom  the 
potential  vortlclty  equation  (1)  using  Identities  (8) 
to  (10)  and  the  expansion  (l6),  we  get  to  first  order. 


.(1)  0'°' 

(20)  ^(f^a,.^,  ,0 


where 


^^^^  dF-  =  ^  ^T+  r  cos  (t)  ^  ^0+  -^-   (  U  +  «   (  ^r 


Carrying  out  the  differentiation  explicitly,  Eq .  (20) 
becomes 

0^°^  ,(1)     0^°^   n) 

(22)      f^  -^  ^  +  f^  (-^)^  w^l^  =  0 
P  P 

where  fp  =  20  cos  ^   and  f^  =  20  sin  (j).   The  first  term 
(v   )  In  Eq.  (22)  Is  the  so-called  g-term.   Next  from 
the  potential  temperature  equation  (6)  using  the  particle 
derivative  operator  (11)  and  the  expansion  (1?)  for 
we  get  to  first  order, 

(25)  0?^  w^^)  =  0  . 


Since  the  atmosphere  at  rest  must  be  stable  we  specify 
that  \u^       >  0;    otherwise  the  perturbation  procedure 


Possibly  this  condition  may  be  relaxed  In  §3-2  (Perturba- 
tion on  Zonal  Flow) . 
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is  invalid.   Then  Eq.  (23)  states  that 

(24)  J-^^    =   0  . 

The  first  order  radial  (vertical)  velocity  is  zero; 

then  Eq.  (22)  states  that  the  (f)- component  of  velocity 

(1) 

V    must  be  zero, 

(25)  v^^^  =  0  . 

Next  from  the  0-component  of  momentum  equation  (12), 
using  the  particle  derivative  operator  (11)  and 
expansion  (16),    we  get  to  first  order, 

(26)  —^ +  f^w^l^  -  f  v^l^  =   0  . 

rp  ^   cos  cp 

But  since  v^-*-^  =  w^"'"^  =  0  from  Eqs .  (24)  and  (25),  this 
geostrophic  relation  degenerates  to 

(27)  -p^^^  =  0  . 

Thus  the  first-order  motion  is  zonal  with  pressure 

p^   ((t),r,T).   Finally  to  first-order,  the  (j)-  and  r-components 

of  the  momentum  equation,  Eqs.  (13)  and  (l4)  respectively 

yield, 

(28)  -4^  +  f^u^i^  =  0 

rp^      -^ 
and 
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r 
P*  '    "       P 

Eq.  (28)  states  that  the  zonal  flow  Is  geostrophic; 


(2?'       -m  (4^'  -  ^^^'  -  v'^'  =  0 


Eq,  (29)  becomes  the  hydrostatic  relation  only  if 
fgU^    =  0,  e.g.  at  the  poles.   Eqs .  (28)  and  (29)  are 
two  linear  equations  for  three  unknown  quantities 
p^  {^,r,r),    p^  i^,r,x)    and  u^  '{^,r,x)    assuming 
p^  ^(r)  and  p^  ^(r)  are  given.   Thus  the  zonal  flow 
is  not  fully  determined  by  the  first-order  equations. 
also  we  note  that  the  first  order  particle  derivative 
operator  (21)  operating  only  on  zeroth  and  first-order 
terms  is  just  the  partial  time  derivative, 
d^-^U    )/dT  =  (  )^,  since  {    )^  =  0   and  v^"^^  -  w^"^^  =  0 
under  these  conditions. 

3)   Second-order  asymptotic  equations  --  Following 
the  same  steps  outlined  in  obtaining  the  first-order 
equations,  the  potential  vorticity  equation  to  second 
order  yields 

p  p 


d 


iil,  lr^(o),(x)..  ^ 


+  dx   ^HToyL^r   ^    +^2"~^  +  ^3^^r  ' 
0(O)p(l) 


p 


■X- 


This  is  exactly  the  same  situation  which  we  find  for  the 
linearized  steady-state  equations,  obtained  without  the 
time  scaling  (cf.  Eckart  [7])- 
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Eq.  (50)  can  be  simplified  since  d^   (  )/dT  =  (  )   and 
the  first  order  relative  vorticity  1!^^    is 


(31)      C    =  1  [vA  '  -  (u'  '  cos  (^)  1] 

r  cos  (|)   w  <f 

(u^  ^  cos  (t))i 
~   ~     r  cos  ^ 

since  v    =0.   The  second-order  potential  temperature 
equation  in  simplified  form  is 


0 


(1) 


(32)  J^^    +TrTrTT  =  0 


0 


Toy 


(2) 

Eq.  (30)  can  be  simplified  by  eliminating  w^    using 
Eq.  (32),  and  the  second-order  potential  vorticity 
equation  becomes 


f^       a  ^  _  ^^'^co3  ^)^      r^eP 


(33)  -^T— -  +  ^  [ -^  + 


r  cos  (p  0J,  'r 


+ 


Finally  the  second-order  momentum  equations  are 

t,4,        f^        +  f  „(2)  .  f   (2)  +  ,(1)  .  0 

rp  ^   cos  9  '^ 

rp^     ^       p^        -^ 

(2)    (u^^h^ 
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The  complete  second-order  asymptotic  equations  are 
Eqs.  (32),  (35), (34),  (3:),  and  (36).   The  first  three 
equations  have  time  derivatives  which  operate  on 
first-order  terms  only.   Hence,  If  the  first-order  zonal 
flow  Is  stationary,  Eqs.  (32),  (33),  and  (3^)  state  that 
the  second-order  motion  Is  also  zonal;   then  the  second- 
order  terms  yield  no  apparently  useful  Information. 
Consequently  In  order  to  obtain  non-trlvlal  second-order 
motions.  I.e.  meridional  planetary  circulations,  we  must 
allow  the  first-order  zonal  flow  to  be  time-dependent. 
Thus  we  see  that  the  firsthand  second-order  asymptotic 
equations  are  strongly  coupled  to  one  another  even 
though  this  coupled  system  of  equations  is  linear.   There 
are  non-linear  first-order  forcing  terms  in  the  second- 
order  equations;  but  from  the  viewpoint  of  the  second- 
order  equations,  the  first-order  terms  are  specified. 
It  is  now  evident  that  we  can  put  very  few,  if  any, 
additional  constraints  on  this  system  of  equations 
without  drastic  effects. 


In  fact,  then  the  motion  is  zonal  to  all  higher  orders 

In  his  careful  scale  analysis  paper.  Burger  [8] 
emphasizes  the  difficulty  in  obtaining  "specific 
InforTTiatlon  on  the  planetary  scale  beyond  expressing 
its  quasi-stationary  character." 
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Now  looking  more  closely  at  the  second-order  equations, 
both  Eqs .  (32)  and  (33)  are  relationships  between  single, 
second-order  quantities  and  time  derivatives  of  first-order 
quantities,  e.g.  Eq.  (32)  gives  (differentiating  with 
respect  to  longitude  0), 

(37)  w^^^  =  0 
and  Eq.  (33)  imposes 

(38)  v^^^  =  0  . 

Eqs.  (37)  and  (38)  Impose  the  following  condition  on  the 
second-order  0-momentum  equation  (3^): 

(39)  p|q^  =  0  . 

(2) 
That  is,  p^  ^(4',t,T)  must  be  Independent  of  0, 

since  Eq.  (39)  precludes  solutions  of  the  form 

p^^^(0,(}),r,T)  =  p^^ho  +  2Tr,(l),r,T)  ;      Eq .  (39)  states 

(2) 

p'    cannot  be  strictly  periodic  in  0 . 

For  completeness  we  give  the  first  and  second  order 
continuity  equations  obtained  from  Eq.  (I5)  and 
expansion  (I6) : 

(to)    pf°'»'^'  +pf°'  dlv  (u(l',vfl',„(l')  =0  . 

Since  v^"^^  =  w^-*"^  =  0  from  Eqs.  (24)  and  (25),  Eq .  (40) 
becomes  Uq   =  0  or  u^  'i^,r,x)    only;  this  gives 
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first-order  zonal  flow  consistent  with  the  other 
first-order  asymptotic  equations.   To  second-order  we 
have 


+p'l'  div  (u(l',vtl',w(l')  =  0  . 


(2) 

This  equation  gives  directly  u^   In  terms  of  first-order 

quantities  and  Is  complementary  to  the  second-order 

(2) 

asymptotic  momentum  equation  for  p^   ,  I.e.  Eq.  (34). 

Thus  the  second-order  continuity  equation  states 
(42)  u^2)  ^  0  _ 

By  Eqs .  (37)  to  (42),  the  second-order  motion  must  be 
Independent  of  0.   To  recapitulate:   the  first-order 
asymptotic  motion  Is  zonal,  geostrophlc,  and  time- 
dependent;  the  second-order  as^nnptotic  motion  is  also 
zonal,  but  in  addition  is  meridional,  non-geostrophlc, 

and  time-dependent. 

(2)  (2) 

Since  now  pA    =  0  (as  well  as  pAq  =0),  the 

second-order  terms  can  be  eliminated  among  Eqs.  (32), 
(33),  and  (34)  and  we  get  a  first-order  time -Invariant, 
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v^eT '_£§_,,  pt°l©w^     .(1) 


+ 


t0^^>rf-^r#T-'r-p'^'il=o 


r       •  -r 

Eq.  (43)  Imposes  an  additional  constraint  on  the 
first-order  solution;  this  condition  (momentum 
conservation)  does  not  seem  to  restrict  the  first-order 
quantities  to  be  steady-state.   However,  the  possible 
physically  interesting  solutions  appear  rather  limited. 
We  remark  that,  although  a  numerical  value  for  the 
time-scale  of  the  fluctuating  first-order  motion  is  not 
determined  from  our  analysis,  the  combined  first-  and 
second-order  system  of  asymptotic  equations  (keeping  in 
mind  the  time-scaling  t  =  et)  leads  us  to  consider  only 
slowly  -  time- varying  motions. 

4)   Initial  and  Boundary  Conditions  —   In  considering 
asymptotic  motions  on  a  planetary  scale,  we  do  not  pose 
ordinary  initial-boundary  conditions  as  we  do  for 
moderately  large-scale  motions  (cf .  Appendix  I) .   Instead 
we  fully  prescribe  the  first-order  time-varying  zonal 
motion,  which  acts  as  a  forcing  function  on  the  second 
order  asymptotic  equations . 
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3-2.      As:>rmptotlc  Equations   Obtained  by  Perturbing   on 

Zonal  Flow,  U{^,r,x) . 

* 

We  can  obtain  a  less-constrained  system  of 

asymptotic  equations  by  perturbing  directly  on  a 
prescribed  zonal  flow,  U((j),r,T);  In  this  case  the 
first-order  asymptotic  equations  yield  time-dependent 
or  steady  state  tesseral  motions.   Thus,  we  can  study 
the  low -wave -number  regime  In  the  neighborhood  of  zonal 
flow. 

Since  the  derivation  of  these  asymptotic  equations 
should  be  clear  from  the  detailed  description  of  the 
algorithm  for  the  simpler,  but  more  constrained, 
asymptotic  equations  given  In  §3-l>  we  merely  summarize 
the  results  here.   The  required  changes  In  the  procedure 
are  the  following: 

1)  the  zeroth  order  solution  In  Eq.  (l6)  Is  replaced  by 

(U,0,0,p^'^\p^°^  0^°h,  which  are  functions  of  ((}),r,T), 

->  -> 

2)  the  zeroth  order  absolute  vortlclty  (20  +  ^  )  Is 

(0,f2-Hn^°^f  +C^°^)  where 

and 

5)    the  particle  derivative  operator  (11)  now  has  a 


^In  the  perturbation  on  the  atmosphere  at  rest,  the  rest 
state  Itself  must  be  prescribed  In  addition  to  the  zonal 
flow. 
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zeroth  order  term, 

d(  )  _     U      MJ_    ,    ^    6^^h     )  ,  ofp2. 
dT    r  cos  (j)  ~^0"  ^  ^   dT    +  ^^^  i 

where  the  first  order  operator  d    (  )/dT  Is  defined 
in  Eq .  ( 21 ) . 

The  zeroth  order  equations,  which  the  prescribed 
zonal  flow  U((j),r,T]  must  satisfy,  are 


(^5)    P^°^  =  0  ,  (  zonal  ) 


U^  tan  (j)    P 


(45)     ^   ^^^  ^   + 


(0) 

^  +  f^U  =  0  ,     ("geostrophlc") 


rp 

(^7)     -  ^  +  -Tot  +^  -  ^2^  =  0  ,    ("hydrostatic") 
p       r 

These  zeroth  order  equations  (45,^6,47)  are  comparable 
to  the  first  order  asymptotic  equations  (27,28,29). 

The  first  order  asymptotic  equations  which  follow 
are  the  principal  approximation  equations.   Prom  the 
potential  vortlclty  equation  (1)  we  get 

(48,    -^^^i—  ^u'-^V^'^0i°'c'^'  1  . 


\^  i 

rp   cos  9 


r    '      ^  r 


.l!.!^    reS^'-^P-,   . 
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+  (fj  +  c'"')  (©r 


(X,   @^°'p'^' 


d 


(1) 


+ 


m 


dT 


P'^L 


(fo-^^^^) 


0 


p 

0) 


(0) 


+ 


r    +tvc'°''0f' 


=  0 


where 


(49) 


m   ^ 
dT 


V 


(1) 


(  )   + 


(  )i  +  ^^^h  ) 


From  the  potential  temperature  equation  we  get 
(50) 


U 


r  cos 


0  . 


From  the  (0^(j),r)  momentiim  equations  respectively  we  obtain 

.(1) 
^(1)   =  0 


(51) 


U  (_^  +  ,(1)  _  v^l^tan  A     -.-^ 

r   Icos  (t)  ^;    V   cos  <t) 


+ 


2        5 


(l)^nv.  A.        1  f     (1)  Pf  'P 


(52)    ^nu^^^tan^^^ 

(1)     1   /(I)   P;-'P 


roT  IP, 

rp^    \^'       p 

(0)J1) 


o)p(i)^ 

roy 


+  f  u^-^^  =  0 


(53) 


2Uu 


+ 


Toy   l^r 


roy 


f^u^^^  =  0 


The  first  order  equations  (48,50,51,52,53)  are  comparable 
to  the  second  order  equations  (33:.  32, 3^,  35, 36)  •   Eq.  (5I) 
st.^tes  that  the  first  order  meridional  motion  is 
non-geostrophic.   Only  equations  (48,50)  have  time 
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derivatives;  and  the  general  structure  of  these  first 
order  asymptotic  equations  are  rather  similar  to  the 
second  order  asymptotic  equations  in  §3.1.   However, 
these  first  order  equations  exhibit  a  coupling  between 
the  meridional  motion  and  the  longitudinal  (0)  variation 
(or  zonal  wave  niimber)  .   In  this  asymptotic  formulation 
the  zeroth  order  zonal  transport  of  the  first  order 
potential  vorticity  and  potential  temperature  is  directly 
related  (proportional)  to  the  first  order  meridional 
transport  of  the  zeroth  order  potential  vorticity 
and  potential  temperature  respectively.   Non-stationary 
planetary-scale  wave  motions  occur  only  with  unsteady 
zonal  motion.   Again  as  in  §3-1  a  thorough  numerical 
study  of  these  asymptotic  equations  will  be  necessary 
in  order  to  test  their  validity  for  describing  either 
steady  or  slowly-varying  planetary-scale  motions  of  our 
atmosphere . 
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Appendix  I.   The  Tangent-Plane  Approximation  for 
Adlabatlc  Flow. 

The  tangent-plane  approximation  Is  easily  obtained 
from  the  derivation  outlined  In  §3-1;  fo^  this  special 
case,  we  get  f^  =  2  0cos  (f)  =   constant  =  fo  and 

^  "=0 

f^  =  20  sin  (t  =  constant  =  f^   .   The  zeroth  order 
terms  are  unaffected. 

Prom  the  potential  vortlclty  equation  for  flrst- 
and  second-order,  respectively,  we  get 

0    p 

Instead  of  Eq .  (22)  and 


(I-l)  f^   (-^)^  w^l)  =  0 


(1-2)   f   ( ^)^  w    +^:p-  i-77TT  L^^    C    + 

op''  P 

.  (01 

p 


0(O)J1) 


-f5(0^"-^^^T^'l'^° 


o  p 

where 

dT   "  ^  ^T  ^  r  cos  (j)  ^  ^0  ^  r    ^  ^^ 

Is  the  horizontal  particle  derivative  operator  since 
w^-'-^  =  0  ;  and 

C^l)  =  ^  [v^^^  -  (u^^^cos  4)).] 

^      r  cos  (|)  '■  9      ^        ^'9' 

Instead  of  Eqs .  (5O),  (21)  and  (3I),  respectively. 
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Prom  the  potential  temperature  equation  for  flrst- 
and  second-order,  respectively,  we  obtain  the  same  equations 
as  for  the  planetary  case,  namely  Eqs .  (25)  and  (32). 
Both  Eqs.  (I-I)  and  (23)  state  that  w^""-^  =  0  for  0^°^  >  0; 

and  from  the  second  order  Eqs.  (1-2)  and  (32)  we  can 

(2) 
eliminate  the  radial  velocity  w^    and  get  a  modified 

(by  the  d©^   /dr  term)  potential  vortlcity  conservation 

equation  of  only  first  order  quantities; 

(7^)      d^^\„(i)^%  r.p^Q^e^^\    (1),   o 

p      ©^ 

The  first  order  momentum  equations  are 

Pq  r-i  ^ 

(1-4)   -j^ -  f   v^^  =  0   (geostrophic) 

p    r  cos  (j)    ^o 

1) 
(0)  .    ^o 


(1-5)   —7^^ +  f-^  u^-"-^  =  0        (geostrophic) 

P 


(1-6)    p^^^  ^^p^^K   0 


We  note  that  the  (f)-component  of  the  coriolis 

parameter  fp  does  not  appear  in  these  first  order 

o 
equations . 

By  introducing  a  new  variable  ^  =  ^ ,  the 

^o 
momentum  stream  function,  and  introducing   cartesian 

coordinates  by  the  transformations  (  )r./^   cos  ^  =    {    )    , 
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(  )^/r  =  (  )y,  (  )r'  =  ^  ^z  ^""^  ga^/r^  %  g  ,  the 

Eqs .  (1-3)  to  (1-6),  inclusive,  become  in  more  natural 

coordinates  for  the  plane  case, 

(l-T)   Vt«+^f(  '!(o)'   'z^^"°° 


(1-8)   p'°'(u'^',  v'-"-')  =  (-'I    ,f    )        (geostrophlc) 

y  X 

(1-9)   P^^^  =  -i^  ^z 


where 


(I-IO)   4?^=  (  ),  +"(!'(  )^+v(l'(  )^ 


(l-ll)       Al-    =    T^   +   ?yy 


,  ,     0'°'       p'°' 

(1-12)   S(°'  =  -^  =  -  (^  ^-  -^)  >  0 


0^"'      p'"'    (o>"')' 


(1-13)   (e<°')2=7RT'°'  =-^:e!^ 

r 

in  which  S    is  the  static  stability  parameter, 
c    is  the  velocity  of  sound,   R  is  the  gas  constant, 
and  T^  ^  is  the  static  temperature;  T^    is  a  given 
function  of  height  z. 
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The  initial  condition  for  the  system  of  equations 
(1-7)  and  (1-8)  is: prescribe  ^(x,y,z,o)  initially. 
The  boundary  conditions  are  pw  =  0  at  z  =  0,oo; 
to  0(e  )  this  condition  becomes  p^   w    =  0  or  by 
Eq.  (32) 

(c^  '  ) 

for  z  =  0 ,  oo  . 

Presiimably  the  geostrophic  potential  vortlcity 
conservation  equation  (1-7)  describes  slowly  varying, 
moderately  large-scale  motions  in  the  atmosphere 
away  from  the  equator.   Again  as  in  the  main  text  of 
this  paper,  numerical  estimates  on  how  slowly  varying 
and  what  horizontal  length  scales  are  involved  is 
missing  from  this  analysis.   Some  information  about 
these  Important  requirements  is  given  in  Appendix  II. 
We  mention  that  we  have  not  been  able  to  include  the 
P-term  (cf.  Eq.(22))  in  Eq .  (1-7)  in  a  consistent  manner. 
To  calculate  the  second  order,  non-geostrophlc  system 
of  linear  equations,  similar  to  that  derived  in  Ref.  [2], 
Is  quite  straightforward. 
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Appendix  II.   Vertical  Structure  of  a  Rectilinear 

Geostrophlc  Vortex  Line  Element  In  Adlabatlc  Flov?; 
Relationship  between  the  Strength,  Horizontal  Length 
Scale,  and  the  Static-Atmospheric  Distribution. 

The  vortex  line  element  Is  obtained  by  seeking 
singular  solutions  of  the  geostrophlc  potential  vortlclty 
conservation  equation  (1-7).   This  concept  Is  strictly 
valid  In  the  tangent-plane  approximation  --  but  also  to 
good  approximation,  we  feel.  In  application  to  flows  In 
middle  latitudes  so  long  as  the  motion  of  these  vortex 
line  elements  does  not  deviate  latltudlnally  very  far  from 
the  Initial  latitude  position.   We  assume  steady-state 
motion  and  seek  solutions  of 

(II-l)     {A?  +-^  [( ^^^^y^ )^  +^^]]   =  5(r)T(z) 

2   2  1/2 
where  &Cr)  Is  the  delta  function  of  r  =  (x  +y  )    ,  I.e. 

5(r)  Is  singular  at  r  =  0  and  Its  Integral  over  the 

entire  x-y  plane  Is  equal  to  a  constant,  say  one; 

r(z)  Is  the  strength  of  the  vortex  line  element,  and  Is 

to  be  determined.   Since  we  assume  that  the  only  flow  In 

the  entire  three-dimensional  space  derives  from  Eq.  (II-l), 

then  a  potential  temperature  invariant  also  exists 


(II-2)  (^z  +  ^J§)^2  ^^T  =  ° 
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We  note  that  In  time-dependent  adlabatlc  flows,  the 
potential  temperature  is  advected  with  the  horizontal 
flow  only  at  the  boundaries,  e.g.  the  ground  (cf. 
Appendix  I).   The  time  derivative  of  Eq .  (II-l),  which 
is  also  a  time  invariant,  combined  with  Eq .  (II-2), 
implies  that  we  are  seeking  solutions  of  a  simpler 
equation  than  Eq .  (Il-l),  namely 

'\ 

(II-3)         (A^  +— 2  1/  )  =  5(r)  .  r(z) 

g    z 

or  using  Eq.  (II-2)  again, 

'I 

(11-^)       (A^ jr^   ^)  =  5(r)   r(z)  . 

Since  Eq .  (II-4)  only  has  horizontal  derivatives  in  the 
X-  and  y-directions  (not  z),  there  is  an  axially-symmetric 
singular  solution  of  Eq.  (II-4)  of  the  form, 


(II-5)  ^0  -  r(z)   Kq[K(z)  .  r] 


where  K(z)  is  given  by 

^3 
(II-6)  K(z)  =  -TTTT^- 

c^°^(z) 

and  K^CKr)  is  the  modified  Bessel  function  of  the 

second  kind.   We  note  that  solution  (II- 5)  is 

not  a  "separation-of -variables"-type  solution  of  Eq.(II-4) 
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A   functional   form  for  r(z)    can  be   found  by 
replacing  Eq .    (II-2)   by 


(II-7) 


hK 


^      + 

2 


17^ 


Y    = 


0 


r(z) 


¥e  note  that  the  right-hand  side  of  Eq .  (II-7)  Is 
regular  in  r.   Then  It  Is  easily  seen  that  r(z)  may 
take  the  form 


(II-8) 


r(z) 


Tq  exp  \-  g 


da 


0  [ctO'(a)f 


where  Tq  =  constant;  and  the  solution  (II-5)  of 
Eq .  (II-4)  also  satisfies  Eq.  (II-7) •   Hence  we  can 
write  solution  (II-5)  more  explicitly  as 

z 


(II-9)      YQ(r,z)    =   TQ-KQEfdzl-r]    exp  l-.&f (off 


2 


where  k(z)  is  given  by  Eq.  (II-6) .   Solution  (II-9)  relates 
both  r(z)  and  k{z)    to  the  given  static  atmosphere;  and 
therefore  there  is  a  functional  relationship  between 
the  strength  r(z)  and  the  horizontal  length  scale 
(<:"  (z),  aside  from  an  arbitrary  constant  Tq  .   From 
Eqs .  (II-6)  and  (II-8)  this  relationship  is 

z 


(11-10) 


r(z) 

^0 


=  exp 


JiE. 


5  0 
^o 


f'^^i^)    da  >  ; 
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or,  differentiating  with  respect  to  z, 

,   dr(z)   d[in  r(z)] 
(^^-11^   r^  ^^  = di =  -  -f-  '^^-^ 

which  eliminates  the  Integration  constant  r^^ .   Since 
the  right  hand  side  of  (ll-ll)  is  definitely  negative, 
dr(z)/d(z)  must  be  negative  when  friction  can  be 
ignored,  that  is,  not  close  to  the  ground  and  possibly 
not  extremely  high  up.   Me  can  get  a  somewhat  simpler 
expression  for  r(z)  than  Eq.  ( 11-10),  using  the 
hydrostatic  condition  which  holds  for  the  static 
atmosphere. 


(11-12)         In  {^^]    =  ^  In  [^--] 


where  7  =  C  /C   and  p^^  is  the  pressure  at  the  ground, 
say.   Eq.  (11-12)  states  that 

(11-13)     r(z)  =  (constant). [p^°hz)]^/^ 

=  (constant)-  p^°^z)  •  0^°^z)  • 

That  is,  the  strength  r(z)  is  proportional  to  the 
product  of  the  static  density  and  static  potential 
temperature. 
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In  addition,  we  note  that  there  is  an  identity  (for 
fixed  z)  which  serves  to  define  the  strength  r(z), 

(11-14)   r(z)  =  ^  ^  ds  -  K^{z)   JJTq  dx  dy 

C  R 

which  might  be  useful  for  calculating  r(z)  from  actual 
data  for  mature  hurricanes  and  typhoons  (cf.  Refs .  [2] 
and  [9]);  the  calculated  results  can  be  compared  to 
Eqs.  (11-10),  (11-11),  or  (II-I3) •   In  Eq .  (II-14) 
C  is  an  arbitrary  horizontal  closed  contour,  surrounding 
the  hurricane  and  enclosing  the  area  R. 
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